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Abstract 

Recent massive numerical simulations have shown that the response of a 
"stochastic resonator" is enhanced as a consequence of spatial coupling. Sim- 
ilar results have been analytically obtained in a reaction-diffusion model, using 
nonequilibrium potential techniques. We now consider a field-dependent diffu- 
sivity and show that the selectivity of the coupling is more efficient for achiev- 
ing stochastic-resonance enhancement than its overall value in the constant- 

diffusivity case. 
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The phenomenon of stochastic resonance (SR) — namely, the enhancement of the out- 
put signal-to-noise ratio (SNR) caused by injection of an optimal amount of noise into a 
periodically driven nonlinear system — stands as one of the most puzzling and promising 
cooperative effects arising from the interplay between deterministic and random dynamics 
in a nonlinear system. The broad range of phenomena — indeed drawn from almost every 
field in scientific endeavor — for which this mechanism can offer an explanation has been put 
in evidence by many reviews and conference proceedings, Ref. [|l|] being the most recent and 
comprehensive one, from which one can scan the state of the art. 

Most phenomena that could possibly be explained by SR occur in extended systems: 
for example, diverse experiments are being carried out to explore the role of SR in sensory 
and other biological functions [0 or in chemical systems 0. Notwithstanding this fact, 
the overwhelming majority of the studies made up to now are based on zero-dimensional 
systems, while most of the features of this phenomenon that are peculiar to the case of 
extended systems — or stochastically resonating media (SRM) — still remain largely to be ex- 
plored. Particularly interesting numerical simulations on arrays of coupled nonlinear oscil- 
lators have been recently reported [^], indicating that the coupling between these stochastic 
resonators enhances the response of the array, which exhibits moreover a higher degree of 
synchronization. This effect has its counterpart in the continuum, as a study on the over- 
damped continuous limit of a (p^ field theory shows ^ . Recently — by exploiting the previous 
knowledge of the nonequilihrium potential (NEP) |^ for a bistable reaction-diffusion (RD) 
model 1^ — one of us has shown analytically that the SNR increases with diffusivity in the 
range explored 0. 

While considering a constant diffusion coefficient D is a standard approach, it is not the 
most general one: it is reasonable to expect that the reported enhancement in the SNR by 
the effect of diffusion could depend in a more detailed way on D. In this regard, see for 
instance 0. In this letter we consider the more realistic case of a field-dependent diffusion 
coefficient D{(f){x,t)), and show that it causes an enhancement of the SNR still larger than 
the one associated with a homogeneous increase of its amplitude. 
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The model under study — a one-dimensional, one-component RD model describing a sys- 
tem that undergoes an electrothermal instability |^ — can be regarded as the continuous limit 
of the coupled system studied by Lindner et al [Q. The field t) might describe the (time- 
dependent) temperature profile in the "hot-spot model" of superconducting microbridges [0. 
This model can be also regarded as a piecewise-linear version of the space-dependent Schlogl 
model for an autocatalytic chemical reaction, and that for the "ballast resistor" , describing 



the so-called "barretter effect" [|10l- As a matter of fact, since in the ballast resistor the 
thermal conductivity is a function of the energy density, the resulting equation for the tem- 
perature field includes a temperature-dependent diffusion coefficient in a natural way ||10||. 



Pointers to other contexts in which a description containing a field-dependent diffusivity 



becomes inescapable have been included in Refs. [pT| , p!2 . 

By adequate rescaling of the field, space-time variables and parameters, we get a dimen- 
sionless time-evolution equation for the field t) 

dt<P{x,t)=d,{D{(P)d,<P) + f{(P) (1) 

where /(0) = —(f> + 9{(j) — 0c) 7 ^'(x) is Heaviside's step function. All the effects of the 
parameters that keep the system away of equilibrium (such as the electric current in the 
electrothermal devices or some external reactant concentration in chemical models) are in- 
cluded in 0c- Moreover, since the value of the field 0(x, t) corresponds in these models to 
the deviations with respect to e.g. a reference temperature Tg > (the temperature of the 
bath) in the ballast resistor or to a reference concentration po in the Schlogl model, it is 
clear that — up to a given strict limit (i.e. = — for the ballast resistor) — some negative 
values of 0(x, t) are allowed. 

As was done for the reaction term a simple choice that retains however the qual- 
itative features of the system is to consider the following dependence of the diffusion term 
on the field variable 

D(0) = Do(l + /l^(0-0e)), (2) 
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For simplicity, here we choose the same threshold 0c for the reaction term and the diffusion 
coefficient. The more general situation is left for a forthcoming work | |13|| . 

We assume the system to be limited to a bounded domain x G [—L, L] with Dirichlet 
boundary conditions at both ends, i.e. 0(±L, t) = 0. The piece wise- linear approximation 
of the reaction term in Eq.(|I]) — which mimicks a cubic polynomial — was chosen in order 
to find analytical expressions for its stationary spatially-symmetric solutions. In addition 
to the trivial solution (j)o{x) = (which is linearly stable and exists for the whole range of 
parameters) we find another linearly stable nonhomogeneous structure (f)s{x) — presenting an 
excited central zone (where (f)s{x) > (pc) for —Xc < x < Xc — and a similar unstable structure 
0u(x), which exhibits a smaller excited central zone. The form of these patterns is analogous 
to what has been obtained in previous related works 0, as is shown in Fig. 1. The difference 
is that in the present case d(f)/dx\x^ is discontinuous and the area of the central zone depends 
on h. 

The indicated patterns are extrema of the NEP, which — among other properties that we 
shall be using presently — is a Lyapunov functional for the deterministic system introduced 
thus far. In fact, the unstable pattern (f)u{x) is a saddle-point of this functional, separating 
the attractors (poi^) and 0s (x) The notion of a NEP has been thoroughly studied, mainly 
by Graham and his collaborators . Loosely speaking, it is an extension to non-equilibrium 
situations of the familiar notion of (equilibrium) thermodynamic potential. For the case of 
a field-dependent diffusion coefficient D{(f){x,t)) as described by Eq. ([l|), it reads |7|,p!0[ 



m = j^' |- D{<P')f{<P') d<P' + \ (^(0)|^) I dx. (3) 

Given that dt(j) = — {\/ D{(j)))6T / 6(j) one finds JF = — / [6T /dcjyf' dx < 0, thus warranting 
the Lyapunov-functional property. This NEP functional offers the possibility of studying 
not just the linear but also the nonlinear — in the case at hand the global — stability of the 
patterns, following its changes as the parameters of the model are varied |^. 

For a given threshold value 0*, both wells corresponding in a representation of the NEP 
to the linearly stable states have the same depth (i.e. both states are equally stable). Figure 



2 shows the dependence of on the parameter 0c- As in previous cases, we analyze only 
the neighborhood of (pc = 0c RW- Here we also consider the neighborhood oi h = 0, where 
the main trends of the effect can be captured. 

Now, with the aim of studying SR, we introduce a weak signal that modulates the po- 
tential around the situation in which the two wells have the same depth. This is accom- 
plished by allowing the parameter 0c to oscillate around 0*: 0c(t) = 0* + <50cCos {Qt + ip), 
with 50c <^ 4>c(t). We also introduce in Eq.(|l]) a fluctuating term ^(x, t), which we model 
(as is customary) as an additive Gaussian white-noise source with zero mean value and a 
correlation function (^(x, t) C,{x', t')) = 27 6{t — t') 6{x — x'), thus yielding a stochastic partial 
differential equation for the random field (f){x, t). The parameter 7 denotes the noise strength 
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As in previous works M , we exploit a generalization to extended systems of the Kramers- 



like result for the evaluation of the decay time or "mean-first-passage time" (r) |1T5|JT6[| . Here, 



those results are extended to the case of field-dependent diffusivity, yielding [13 



f W[0,0.] ] 

(r) = roexp| — — — |. (4) 

The functional W[0, 0j] (0j indicates the initial metastable state, which at each instant may 
be either 0o or 0^), that is the solution of a Hamilton- Jacobi-like equation, has the following 
expression 

with f/(0) = Jq d(f)' f{(j)'). The prefactor tq in Eq.(^ is essentially determined by the curva- 
ture of the NEP JF[0, 0c] at its extrema. 

The calculation of the SNR proceeds, for the spatially extended problem, through the 
evaluation of the space-time correlation function {(f){y,t)(j){y',t')). To do that we use a 
simplified point of view, based on the two-state approach [0, which allows us to apply 
some known results almost directly. To proceed with the calculation of the correlation 
function, we need to evaluate the transition probabilities W± oc (t)^^, which appear in the 
associated master equation. For small 50c, 



Solving such a master equation up to first order in 5(f)c it is possible to evaluate the correlation 
function. Its double Fourier transform, the generalized susceptibility S{k.,uj), factorizes 
in this approach, and the relevant term becomes a function of uo only (the corresponding 
expressions are omitted, see [0 for details). 



It is worth noting that many of the results exposed here (e.g. the profiles of the station- 
ary patterns and the corresponding values of the NEP) are exact. The only approximations 
involved in the calculation of the SNR are the standard ones, namely the Kramers-like ex- 
pression in Eq.(^) and the two- level approximation used for the evaluation of the correlation 
function . 



Using the definition from Ref . [|T^ for the SNR at the excitation frequency (here indicated 
by -R), the final result is 

i?^(A)2exp(-n;[0,0.]^./7), (6) 

where A = ((i>V[0, 0i]/(i0c)0j'^0c, and tq is given by the asymptotically dominant linear 
stabihty eigenvalues: tq = 27r(|A""| A"**)^^/^ (A™ is the unstable eigenvalue around 0^ and 
A"** is the average of the smallest eigenvalues around 0o and (pg)- Equation (^ is analogous 
to the results in zero- dimensional systems, but here A, tq and W[0, ^i],/)* contain all the 
information regarding the spatially extended character of the system. 

In Fig. 3 we depict the dependence of R on the noise intensity 7, for several (positive) 
values of h. These curves show the typical maximum that has become the fingerprint of the 
stochastic resonance phenomenon. Figure 4 is a plot of the value Rmax of these maxima as 
a function of h. The dramatic increase of Rmax-, of several dB for a sma// positive variation 
of h, is apparent and shows the strong effect that the selective coupling (or field- dependent 
diffusivity) has on the response of the system. 

The present prediction prompts to devise experiments (for instance, through electronic 
setups) as well as numerical simulations taking into account the indicated selective coupling. 



This result could be of relevance for technological applications such as signal detection and 
image recognition, as well as for the solution of some puzzles in biology (mammalian sensory 
systems, ionic channels in cells). 

The present form of analysis is being extended to (the bistable regime of) multicomponent 
models of the activator-inhibitor type since — in addition to their applications to systems of 
chemical (e.g. Bonhoffer-Van der Pol model) and biological (e.g. FitzHugh-Nagumo model) 
origins — these models are related to spatio-temporal synchronization problems An 
effective treatment of models of this type gives rise to a non-local coupling, which would 
compete with the nearest neighbour coupling D{(j)) presented here ||18|| . 
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FIGURES 

FIG. 1. (f)s{y) is the stable pattern, and (puiv) the unstable one (similar in form but exhibiting 
a smaller excited central zone); they are extrema of the NEP: 4>uiy) is a saddle-point of this 
functional, separating the attractors (f)Q{y) and (f)s{y) [7]. Here we indicate (f)o{y) and (f)s{y) with a 
full line, ^^(y) with a dotted line, while the value of 0c is shown with a dashed line. 

FIG. 2. Dependence on 4>c of the value of F calculated for the patterns in Fig.l. The stable 
branch Fg (indicated with a full line) and the unstable one F„ (indicated with a dotted line) collapse 
at a critical value of 0c- At a certain value (pc = (pc (indicated by an arrow), the value of Fs becomes 
positive and 0s (x) becomes metastable. 

FIG. 3. SNR i? as a function of the noise intensity 7 (Eq.(6)), for three values of h: h = 0.0 
(full hne), h = —0.25 (dashed line) and h = 0.25 (dotted line). We have fixed L = 1, Dq = 1, 
S(pc = 0.01 and n = 0.01. 

FIG. 4. Maximum Umax of the SNR curve (Fig.3) as a function of h, for three values of Dq: 
Dq = 0.9 (dashed line), Dq = 1. (full line) and Dq = 1.1 (dotted line). The arrows a and b 
indicate the response gain due to an homogeneous increase of the coupling and to a selective one 
respectively. The larger gain in the second case is apparent. The inset shows the dependence of 
Rmax on Do for h = —0.25 (lower line), ^ = and h = 0.25 (upper line). 
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